On the mechanism of the highly viscous flow 
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The asymmetry model for the highly viscous flow postulates thermally activated jumps from a 
practically undistorted ground state to strongly distorted, but stable structures, with a pronounced 
Eshelby backstress from the distorted surroundings. The viscosity is ascribed to those stable dis- 
torted structures which do not jump back, but relax by the relaxation of the surrounding viscoelastic 
matrix. It is shown that this mechanism implies a description in terms of the shear compliance, 
with a viscosity which can be calculated from the cutofi' of the retardation spectrum. Consistency 
requires that this cutoff lies close to the Majcwell time. The improved asymmetry model compares 
well with experiment. 

PACS numbers: 64.70.Pf, 77.22. Gm 
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I. INTRODUCTION 

There are many theories and models of the flow 
processi"— , but two central points are still not clear. The 
first is the detailed development of the viscosity from the 
thermally activated jumps between inherent states^ii^^ 
in the highly viscous liquid. The second is the role of the 
elastic interaction between different jumps. 

The phenomcnological asymmetry modeP"^ focuses 
precisely on these two central points. It is able to ex- 
plain the asymmetry observed in aging experiments at 
the secondary relaxation peak^. In these experiments, 
one gets a sudden small rise of the peak height after a 
small temperature rise even before the energy landscape 
adapts itself to the new temperature. This rise is only 
understandable if there is an average energy difference 
of about 4 ksT between the structural minima of the 
secondary relaxation. 

The asymmetry model ascribes the asymmetry to a 
structural rearrangement of an inner core (a "coopera- 
tively rearranging region"^), consisting of N atoms or 
molecules, with an N of the order of hundred. The struc- 
tural rearrangement leads to a change of shape, which in 
turn leads to an elastic distortion, both of the core and 
of the surrounding matrix, an elasticity theory problem 
solved by Eshelby^. The Eshelby concept, first intro- 
duced into the field to explain the plastic deformation of 
metallic glasses^and closely related to the concept^" of 
shear transformation zones, is a good starting point for 
a quantitative description of the flow process, not only 
because it provides a quantitative relation between the 
local relaxational jump and its coupling constant to an 
external shear, but also because it allows for a simple 
mean-field treatment of the interaction between different 
relaxing entities. 

According to the asymmetry model, the flow proceeds 
by a jump from the unstrained ground state over a rel- 
atively high barrier into a distorted state. If the barrier 
is high, the inner domain remains for a long time in the 
distorted state. If this time is long enough, the inner core 
can relax by relaxation processes in the neighborhood, to 



become a new unstrained ground state. 

If this is in fact the appropriate physical basis, the cal- 
culation of the shear response in the former treatment^ 
can be considerably improved by relating the viscosity to 
the cutoff of the retardation spectrum. This is the aim 
of the present paper. 

The following section II begins with a reformulation of 
the asymmetry model in terms of the shear compliance, 
deriving expressions for the viscosity and the retardation 
spectrum. Section III describes fits to experimental data 
from the literature. Section IV discusses the success and 
the limitations of the treatment, section V contains the 
summary. 



II. THE COMPLIANCE FORMULATION OF 
THE ASYMMETRY MODEL 

One has two equivalent textbook descriptions of the 
shear spectrum of a liquid^^, a relaxation spectrum H(t) 
for the description of the complex shear modulus G(w) 
(t relaxation time) 



H{T)iujT 
1 + iujT 



d Inr 



(1) 



and a retardation spectrum L{t) for the description of 
the complex shear compliance 
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Here Goo is the infinite frequency shear modulus and 
T] is the viscosity. A third material constant hidden in 
this equation is the recoverable compliance , the elas- 
tic compliance plus the integral over the retardation pro- 
cesses 
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Goo 
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(3) 



Eq. ^ makes a separation of two independent con- 
tributions to the compliance, the retardation spectrum 
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and the viscosity. The retardation spectrum is due to 
back-jumps into the initial inherent state, the viscosity is 
due to no-return processes. The retardation description 
separates these two influences, the relaxation description 
does not. In the relaxation description of eq. (1), 77 is 
given by 



V = 



TH{T)d Inr 



(4) 



and is an inherent property of the spectrum. 

In principle, one could also separate H{t) into a flow 
and a relaxation component, but the truncation would 
be far less obvious than in the L{t) case. A formally 
similar problem is encountered in the field of dielectric 
relaxation, where the dielectric modulus could also in 
principle include a zero frequency component. But the 
usual way is to add an independent conductivity term to 
the dielectric susceptibility. 

Since the asymmetry model^^ predicts an independent 
flow process (though mediated by the retardation pro- 
cesses), it should be formulated in terms of the com- 
pliance description. This was not done in the original 
derivation^"^ , but will be done here. As a further im- 
provement, the Poisson ratio ^'^'^^ is explicitly taken into 
account. 

The asymmetry modeP^ considers an embedded core 
volume of N atoms or molecules (atomic or molecular vol- 
ume Va) which can undergo structural rearrangements. 
The shape difference of the inner volume Nva is either a 
shear (the shear angle between the two stable struc- 
tures of the inner volume in radian) or a volume change 
dv or a mixture of both. For a shear and an infinitely 
hard surrounding medium, the asymmetry (the energy 
difference between final and initial state) A would be 
GooNvae'^ /2, where Goo is the (infinite frequency) shear 
modulus. But we consider a surrounding medium with 
the same elastic constants as the inner part. In this case, 
the medium distorts as well and one gets a lowering of 
the elastic energy by about a factor of two, slightly de- 
pending on the Poisson ratio ly of the substance^^ 
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where for a spherical inner volume the Eshelby-Poisson 
reduction factor jp is given by 



7p 



7-5v 
15(1 - v) 



(6) 



and the coupling to an external shear e in the direction 
of Ci is given by 



dA 

-g^ = IpGooNVaCi 



(7) 



This implies that the influence of the Poisson ratio can be 
taken into account by replacing GooVa in the equations 
of the earlier paper^^^ by l^pGoo^a- 



A given asymmetry A defines the surface of a sphere 
in the six-dimensional distortion space. For a constant 
density of stable structures in distortion space, one gets 
a number of states n(A)c?A between A and A -f dA 



n(A) 



cn 



m (27pGoo«a)3' 



(8) 



where cat is a dimensionless and temperature- 
independent number. The partition function 
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Cn 



(27pGooiV«a)3 



A^e'-B^dA 
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is easily evaluated. With the definition of the character- 
istic multi-minimum parameter /at 



J. Cn 



the partition function reads 

Z=l + 2fN, 



(10) 



(11) 



where the 1 stands for the ground state and the 2/jv for 
the distorted excited states. 

The second order term results from the integral 



fNc'^lpG^NVa 



where the function P5g(A) is given by 
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Since this function is the same as in the treatment 
without the Eshelby-Poisson factor 7p, the conclusion of 
an average asymmetry of -iksT remains unchanged. 

To second order, the partition function with an applied 
strain e is 



Z -1 + 2/w H j-^ e 

kuT 



(14) 



With eq. ([T4|) . one can calculate the contribution 
5Gn of the volume Nva to the reduction of the shear 
modulus, using the free energy F — — ksT In Z and 
SGn = {l/Nva)d^F/de'^ 
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Go 
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This is the reduction of the shear modulus without inter- 
action between different structural rearrangements. 

Let us next consider the infiuence of the dynamics, 
characterized by the variable v — IuTi,, where r„ is the 
relaxation time. Let us denote by foiy) the barrier den- 
sity of structural rearrangements of the N molecules or 
atoms, in the sense that —SGn/Goc from these relax- 
ations (without interaction) between v and v + dv is given 
by fo{v)dv. 
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With this definition, it turns out to be surprisingly 
easy to take the elastic interaction between different re- 
laxation centers into account. At the time r^,, the shear 
modulus is reduced by ah the relaxation centers with 
shorter relaxation times than to the value G^. For 
a given relaxation, this means that its asymmetry A is 
reduced from its initial value by the factor Gy/Goo (as- 
suming that the inner domain stress energy decreases by 
the same amount as the stress energy outside) . Since the 
total effect for all relaxations between v and v + dv is 
an integral over the full distortion space, one gets again 
the integral over the function psai^) of eq. (fT^ . If all 
A-values decrease by Gy/Goo, the integral increases by 

(Goo/Gy) . 

On the other hand, the coupling constant to an exter- 
nal strain e also decreases by the factor Gy/Goo- Since 
the coupling constant enters with its square into the de- 
crease of the shear modulus, the barrier density f{V) 
with interaction must be related to the barrier density 
fo{V) wit/ioitt interaction by 



(16) 



The solution to this equation has been given earlier^®. 
For a given fo{v), f{v) reads 



fiv) 



foiv) 



2/3 ' 



(17) 



a solution which brings the shear modulus down to zero 
if fo(v) integrates to 1/3 (the 1/3- rule). 

But the essential point of the present paper, the one 
that was not understood in the previous attempt o^^'^^ , 
is that the shear modulus does not go down to zero; it 
goes down to 1/ Jg and stays there, while the independent 
viscosity processes do the rest. 

This being so, one needs to formulate the problem in 
terms of the retardation spectrum L{t) and the viscosity 
rj. Let us do that. 

Since the relaxation centers do not couple to the stress, 
but to the strain (see eq. ([7])), one has to formulate 
the description of their effect in terms of a relaxation 
spectrum like eq. ([1]), but with a finite rest modulus 
1/Jg at the frequency zero 



G.(w) = jo 



Go 



f{v)iujTy 

g 1 + iUJTy 



dv. 



(18) 



Gr (w) has to be inverted to obtain the shear retardation 
spectrum. Equation (fT8)) implies that the integral over 
f{v) is limited to 1 — l/GooJe- Translating this into a 
limitation of the integral over fo{v) 



fo{v)dv 
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(GJ0)3 



(19) 



Obviously, the function fo{v) must have a cutoff at a 
critical u^, with a critical relaxation time Inr^ = Vr. 



Denote the inherent relaxation density (the one with- 
out cutoff) by fi{v). Let these inherent relaxations decay 
with the time constant Tc. Then the cutoff function for 
fi{v) is given by Tc/ijy -I- Tc). The function /o('u) must 
result from the multiplication of the inherent relaxation 
density ]i(v) with this cutoff function squared 



/o(^') 



(20) 



because only a fraction T^/iry -I- Tc) of fi{v)Tc/ (Ty + Tc) 
returns to contribute to the retardation. The rest, a frac- 
tion Ty/ {Ty + Tc) of fi{v)Tc/ {Ty + Tc), decays via the vis- 
cosity escape. Its contribution is only sizable if Ty ~ Tc- 
This opens up a way to calculate the viscosity. 

For simplicity, the asymmetry energy A as well as all 
other energies will be counted in units of ksT in the 
following considerations. 

To get the viscosity, assume a stationary fiow with an 
applied shear stress a. The "back-lag" e?, of the equi- 
librium strain is crJg. Take a relaxation center with an 
asymmetry A which happens to have its inner structural 
shear displacement in the direction of the applied shear 
stress. In the stationary fiow, the population p of the dis- 
torted state is increased by the factor 
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dA 



l + aJ^^iA^pGooNva, (21) 



where equs. ([5]) and ([7]) have been used. 

If the surroundings of the distorted state with energy 
A relax, it becomes a new unstrained ground state. In 
the transition, the volume Nva flows by the shear angle 
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IpGooNVa 



One then gets a contribution to the flow 



5e = 2Acr j; 



(22) 



(23) 



Note that the number TV, the Eshelby-Poisson factor and 
Goo disappear in the product. If one has to replace Goo 
by Gy, the equation remains valid. 

One gets the same contribution from a relaxation cen- 
ter oriented in the opposite direction, because both the 
population change and the flow direction change sign. 
For the generally oriented defect, one gets a factor of 1/6 
from the six components of the distortion, one compres- 
sion and flve shears. 

Integrating over A, one gets a factor of 3 (a factor of 
6 from the integral over pgQ and a factor of 1/2 from the 
normalization). Summing up over all relaxation centers 



IpTcjQ (Tc+Tu)2 



!t{v)dv, 



(24) 



where the 7p in the prefactor results from the weight of 
the excited states in the partition function, equs. (jl4p 
and p^ . 
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Since e = cr/77, this supplies an equation for the 
Maxweh time tm ~ v/Goo 
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TM 



IpTc 



'-fa{v)dv. 



(25) 



The same integral supplies the lifetime Terr of the co- 
operatively rearranging region, because a fraction r^/r^ 
of fo{v) decays with the rate 1/tc. Therefore 



GooJe^M ■ 



(26) 



This lifetime is considerably longer than Tc, because only 
the excited states decay and even there, the decay com- 
petes with the back-jump probability. 

The weight factor fgs of the undistorted ground state 
in the canonical ensemble can be calculated from equs. 
([TT|) and ([13. One finds 
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1 - l/{G^J^,f 
37p 



(27) 



a value close to 1/3. 

In order to achieve consistency, the decay rate I/tc 
should be ascribed to the viscous flow outside the inner 
core. As pointed out by Eshelby^^, one has a fraction 7^ 
of the total elastic energy in the core and a fraction 1 — 7p 
outside. This implies a decay rate for the asymmetry A 
by the viscosity outside 



A 
A 



1 - 7p 



(28) 



The distortion lifts the energy of the distorted state by 
A. The energy increases quadratically with the distance 
from the ground state in distortion space. Therefore it 
lifts the energy at the barrier between the two states by 
A/4, and the relative change of Ty is given by 



3A 1 



7p 



Tm 



(29) 



Tc was defined as the relaxation time at which an equal 
number of inherent states goes into both channels, the 
viscosity and the retardation. This implies the condition 
7\, = 1 at the decay time constant Tc, which should be 
obeyed at the average value A = 4, so the consistency 
condition reads 



TM 



= 3(l-7p). 



(30) 



This means Tm is factor of 3/2 larger than Tq if the 
Eshelby-Poisson factor is exactly 1/2. In any case, Tc 
lies close to the Maxwell time tm- 



III. COMPARISON TO EXPERIMENT 



frequency shear modulus Goo, the viscosity 77, the re- 
coverable shear compliance , the Eshelby-Poisson fac- 
tor 7p and a stretching parameter 7 (the analog of the 
Kohlrausch /3). 

For a problem-adapted fitting procedure, one chooses 
Goo, 7p, Goo^e and Vc ~ In Tc as four of the five indepen- 
dent parameters. Vc stands for the viscosity, because one 
can calculate t] via the Maxwell time from eq. ([50]) . 

The stretching parameter 7 is introduced into the 
equation for fo{v) 



fo{v) = aexp(7(t; - Vc) + j2iv - Wc)^) 



(Tc+T„)2 



(31) 



where the factor a and the curvature 72 are not free pa- 
rameters, but fixed by equs. (fT9|) and (|25|) . If 72 is nega- 
tive, the function is a gaussian with a cutoff at Vc = In Tc. 

If one has a pronounced secondary peak, one needs 
the position and the width of this peak as additional 
parameters. In this case, the height of the peak takes 
the role of the curvature parameter 72, which is set to 
zero. 

The most stringent test for the model is supplied by the 
compliance measurements of Plazek and coworkers^S^— , 
because these measurements separate retardation and 
viscosity explicitly. Fig. 1 shows the fit of the recoverable 
compliance data (the recoverable compliance J{t) — t/rj is 
only due to the retardation) for five molecular glass form- 
ers and for boron trioxyde. In boron trioxyde, the shear 
modulus Goo of the fit was taken from Brillouin"^'^ and 
ultrasonic data'^'*. As it turns out, the recoverable com- 
pliance data do not fix the Eshelby-Poisson-factor 7p very 
accurately, so one can still choose it rather freely. The 
choice is nailed down by the measured viscosity, which to- 
gether with the recoverable compliance fixes the Eshelby- 
Poisson parameter with an accuracy of about 1 to 2 %. 
This joint fit requires a non-zero curvature 72, which is 
sometimes positive and sometimes negative. 

Looking at data from the same substance at different 
temperatures, one can check whether the fitted Eshelby- 
Poisson parameter is temperature dependent. In the case 




1 2 

log(t/s) 



Even with the equations derived in the preceding sec- 
tion, one is left with five unavoidable independent param- 
eters with a well-defined physical meaning: the infinite 



FIG. 1: Asymmetry model fits of the recoverable compliance 
of six glass formers at the glass transition^Si"^. Notation see 
Table I (ARO is Aroclor # 1248 and OTP is orthoterphenyl). 
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of Aroclor^S, in Fig. 2, one finds no visible temperature 
dependence even within the small experimental error. 

The error bars are slightly larger (about 5 to 10 %), if 
one fits measurements of G{uj). Fig. 3 (a) shows propy- 
lene glycol^ and Fig. 3 (b) tripropylene glycol^^ as ex- 
amples for substances without and with visible secondary 
relaxation peak. In the case of tripropylene glycol one can 
take position and width from dielectric measurements'^^ 
and replace the parameter 72 by the height of the peak, 
so one has again five parameters. 

The fit of propylene glycol in Fig. 3 (a) provides a neg- 
ative 72 — —0.012. Setting this value constant at all mea- 
sured temperatures, one gets a practically temperature- 
independent 7 — 0.10. This implies that one sees a gaus- 
sian in fo{v) with its maximum at + 4.1 at all temper- 
atures. In other words, the maximum of the barrier den- 
sity fo{v) shifts with temperature in such a way that it 
keeps the constant relaxation time ratio Tmax/Tc = 64.5 
at all temperatures, where Tmax is the relaxation time 
at the maximum. The same is found for the data^^ 
of triphenylethylene, only with the slightly larger ratio 

'^max I — 181. 

This kind of peak shift is expected in the elastic 
models^*^, which postulate a proportionality between 
the barrier heights and the infinite frequency shear 
modulus Goo- What does not fit is the temperature- 
independent curvature 72, which should increase with 
T^/G^ with increasing temperature. Such a decrease 
can be excluded within experimental error in both cases; 
the fits rather indicate a broadening with increasing tem- 
perature. The finding is relevant for the question whether 
the temperature dependence of Goo can indeed account 
for the full fragility'^^. 

The model gives excellent fits. In fact, the propylene 
glycol fit in Fig. 3 (a) is by no means the best, but 
belongs to the poorer ones of all those done to obtain the 
data in Table I. 

Table I compiles the fitted values of GooJe s-nd the 
Eshelby-Poisson factor 7p. In those cases where mea- 
surements at different temperatures were available, the 
Aroclor result of a temperature-independent 7p was con- 
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firmed within experimental accuracy. The recoverable 
compliance increased with increasing temperature in 
the two cases dibutylphthalate and tripropylene glycol, 
both substances with a secondary relaxation peak, and 
decreased with increasing temperature in glycerol and 
meta-toluidine. 

The values of 7^ of Table I look very reasonable. 
The value 0.57 found in orthoterphenyl is slightly larger 
than the one calculated from eq. ([6]). The longitudi- 
nal and transverse sound velocities of orthoterphenyl at 



its glass transition temperature are 



43 



vi 



2.94 km/s 



and vt = 1.37 km/s, respectively. The Poisson ratio 
= {{vi/vt? - 2)(^/wt)' - l)/2 is 0.361, so 7p = 0.542. 
One expects (and finds) a similar value for all the other 
molecular glass formers in Table I. 

The value 0.48 of silica fits its low Poisson ratio of 0.12, 
while the value 0.88 for B2O3 is clearly too large to be 
explained by any possible Poisson ratio in eq. ©. But 
there is a loophole for escape: The Eshelby-Poisson factor 
depends also on the shape of the core region. The point is 
discussed extensively by Eshelby^^. To find such a large 
7p as in B2O3, one needs two parallel extended planes 
shearing against each other in the core region. This could 
indeed happen in B2O3, which tends to form covalent 
planes. Some of the comparatively large 7p-value of OTP 
and other aromatic compounds in Table I could be due 
to the same reason. In the same context, the value 0.5 for 
the two metallic glasses is interesting (from their Poisson 
ratio, they resemble the molecular glass formers). 0.5 is 
the exact valuers, for a needle lying either along x or y in 
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FIG. 2: Asymmetry model fits of the recoverable compliance 
of Aroclor # 1248 at different temperatures^ii. 



FIG. 3: Fits of G{u}) in terms of the new version of the asym- 
metry model in a) propylene glycol b) tripropylene glycol. 
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77". 

silica^^ 


0.48±0.1 


2.0 


boron trioxyde^ 


0.88±0.01 


3.6 


ZresAlT.sCuiT.sNiio^ 


0.50±0.03 


2.5 


Pd4oNi4oP2(7ia 


0.50±0.04 


3.6 


orthoterphenypffi 


0.57±0.01 


2.61 


Aroclor # 12482^ 


0.59±0.01 


2.36 




0.55±0.01 
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TNB-2^ 


0.56±0.01 


3.5 




0.60±0.01 
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ppg36 


0.59±0.05 


2.35 


glycerol^i^ 


0.57±0.05 


6.5-3.5 


propylene carbonate^ 


0.51±0.1 


3.1 


DGEBA^ 


0.56±0.05 


3.0 


m-toluidine^ 


0.54±0.05 
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propylene glycol^ 


0.56±0.05 


5.7 


diethyl phthalate^ 


0.54±0.05 


3.1 


dibuthyl phthalate^ 


0.54±0.05 
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0.55±0.05 


2.5 


triphenylethylene^ 


0.58±0.05 


2.6 


tripropylene glycol^ 


0.52±0.05 


3.1-4.3 



TABLE I: Fitted Eshelby-Poisson factors and recover- 
able compliances for twenty glass formers. Aroclor # 
1248 is a chlorinated biphenyl of the company Monsanto; 
TNB-1 is l,3-bis(l-naphtyl)-5-(2-naphtyl)benzene; TNB-2 is 
l,3,5-tris(l-naphtyl)benzene; 6-PE is bis(m-(m-phenoxy phe- 
noxy)phenyl)ether, a chain of six phenyl rings joined by oxy- 
gens; PPE, polyphenylene, is a commercial diffusion pump oil 
which is chemically very close to 6-PE; DGEBA is an epoxy 
resin, diglycidyl ether of bisphenol-A; DC704 is a silicon oil 
for diffusion pumps. The references denote the sources of the 
experimental data. 



an ea;y-shear (in this case, the needle is neither shortened 
nor lengthened by the shear) . The finding agrees with the 
suggestion^ of a flow in metallic glasses due to the strings 
seen first in numerical work on relaxation in glasses^ and 
later in numerical work on the flow process in liquids^. 



IV. DISCUSSION 

A first question is whether one can indeed do statistics 
with the small number of inherent states of a " coopera- 
tively rearranging region"—, a relatively small volume of 
about hundred atoms or molecules. A recently developed 
four-point correlation methodic allows to estimate the 
size of this region from the temperature dependence of, 
say, the dielectric susceptibility "'^^^^ . The average num- 
ber of inherent states in such a cooperatively rearranging 
region can be calculated from the asymmetry modeP-, 
because the for just two inherent states with asym- 
metry 4k bT in the core is 0.012. To come to the mea- 
sured values of GaoJe of 2 to 3, one needs Jn ~ 0.7. .0.9, 
so there are sixty to eighty inherent states in a single 
cooperatively rearranging region. The size of the region 
increases strongly with decreasing temperature^^—, but 



the number of inherent states remains fixed by the recov- 
erable compliance GooJe^ which according to the results 
in Table I tends to stay constant. The dramatic increase 
of the volume needed for a constant number of inher- 
ent states with decreasing temperature is consistent with 
the well-known dramatic entropy loss which leads to the 
extrapolated Kauzmann catastrophe. The ensemble of 
sixty to eighty states, though small, is large enough to 
make statistical considerations and to calculate a canon- 
ical partition function. 

A further possible objection is that there is no such 
clear separation between outside and inside as assumed 
by the model; the next jump can occur partly inside and 
partly outside of the core region of the former jump. But 
the definition of the core region requires the unstrained 
ground state, which bundles the sixty to eighty jumps to- 
gether (this is also the answer to another possible ques- 
tion, namely why there is no broad distribution of re- 
values). 

The traditional theoretical procedure to calculate 
the shear response^siS, makes use of the Zwanzig-Mori 
formalisn*^, considering appropriate projection opera- 
tors which project the interesting variable out of the full 
set of bath variables and calculating their development 
in time via the Liouville operator. The appropriate vari- 
able in the shear case is the shear strain. Therefore the 
bath should be described in terms of a sum of relaxators 
reacting to the applied shear strain, equivalent to the 
f{v) of the present paper (the Eshelby concepli^S, leads 
to the same conclusion). The viscosity and the recover- 
able compliance are obtained from the first and second 
moments of this relaxation spectrum in relaxation time. 

The analysis of the present work reveals the danger of 
the traditional approach. It induces one to believe that 
the relaxation density which describes the shear response 
reflects the bath processes via their coupling to the shear 
strain. According to the present analysis, this is only 
true in the initial stage, but is the wrong interpretation 
of the final loss of reversibility in the viscous flow. The 
final loss of reversibility requires the detailed analysis of 
the point in time where the system decides to keep or to 
leave a local structure. As a consequence, the viscosity 
is an independent process. This forces one to go over to 
the compliance description of eq. ([2]). 

The present work shows that the Maxwell time marks 
the boundary between relaxation processes owing their 
long relaxation times to a high barrier (on the fast side) 
and relaxation processes owing their long relaxation time 
to the large number of relaxation processes which they 
require (the Maxwell time is itself the first example of 
this second class). At the same time, it is the boundary 
between back-and-forth motion and irreversible flow (the 
paper by Zwanzig^'^ aims at a correct treatment of exactly 
this crossover). 

There is a sizable number of recent numerical at- 
tempts to flnd out what happens microscopically on the 
viscosity time scaled"—. On the accessible relaxation 
time scale of nanoseconds, one finds a dynamic correla- 
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tion length via a four-point correlation functionJ^, which 
seems to diverge with diverging relaxation time accord- 
ing to the Adam-Gibbs recipe^. The present work iden- 
tifies this correlation length with the core diameter of 
an Eshelby rearrangement^'''. Studies of the shear stress 
correlationS2j53 show the strongly nonlocal character of 
the viscous process, in good agreement with the cen- 
tral postulate of the present paper. It is not quite clear 
whether the diverging correlation length dominates the 
viscosity^ or whether the viscous flow requires even more 
spaced, but in any case faster relaxations seem to occur 
on smaller length scales-^. 

The advantage of the numerical work is that one can 
look at details of the atomic motion which are not accessi- 
ble to experiment. The work of the swiss group^ focuses 
on the correlation length for particles with a parallel dis- 
placement vector in deformation-induced transitions be- 
tween inherent states. Consider a pure shear Eshelby 
rearrangement where the core expands in the x-direction 
and contracts along y. The displacements are four mov- 
ing clouds of particles, each of them partly inside and 
partly outside of the core, two moving outside along x 
and two moving inside along y, with a diameter similar 
to the one of the core itself. At the same time, there is 
a density excess in the x-clouds and a density diminu- 
tion in the y-clouds, visible in the four-point correlation 
of the density. Thus one would expect the displacement 
correlation length^ to be close to the dynamic density 



correlation one^, as one indeed finds. 



V. SUMMARY 

The present work describes an attempt to analyze the 
microscopic processes at and around the Maxwell time. 
The analysis bases on the asymmetry model2^, which as- 
sumes relaxational jumps from one unstrained ground 
state to many strongly distorted states of an inner core 
of about hundred atoms or molecules (a "cooperatively 
rearranging region"—) . 

The analysis relates the viscosity to the long-time cut- 
off of the retardation processes. The viscosity is ascribed 
to distorted states staying for such a long time that 
the surroundings adapt, making them new undistorted 
ground states. 

In this picture, the Maxwell time rj/Goo {rj viscos- 
ity, Goo infinite frequency shear modulus) is the border- 
line between reversible faster processes, with a relaxation 
time determined by a thermally activated jump over a 
barrier, and irreversible slower processes, which attain 
their long relaxation times by the necessity to jump over 
many barriers. The viscous flow itself belongs to this 
second category, which forces one to use a compliance 
description. "With compliance comes comprehension" - 
Don Plazek was right'^^. 
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